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1. Introduction 

Heavy ion physics opens the unique opportunity to explore the nuclear 
equation of state (EOS) far away from saturation, i.e. at high densities 
and at non zero temperatures. Investigations of nuclear matter under such 
extreme conditions are crucial for the understanding of the universe, its 
evolution and the formation of the elements, and the evolution of massive 
stars, supernovae and neutron stars. 

The intensive study of heavy ion collisions has became possible by the 
introduction of new generations of accelerators and detectors, which allow 
to study collisions of heavy nuclei from Fermi energies up to high relativis- 
tic energies of a few hundred GeV per nucleon with an almost complete 
characterization of the reaction products of the space-time events. 

Inspite of the theoretical efforts over the last three decades (see ref. [1]) 
the determination of the EOS in heavy ion reactions is still an object of 
current debate. Comparisons with experiments seem to favor a nuclear EOS 
with an incompressibility of K ~ 230 MeV and with a moderate density 
dependence at high densities (soft EOS). In particular, a difficultly in the 
determination of the nuclear EOS is the fact that the nuclear matter is to a 
large extent far away from local equilibrium in a heavy ion collision. These 
non-equilibrium effects are a main feature of energetic nuclear collisions, 
and govern the dynamics of the reaction with relaxation times comparable 
to the compression phase of the process [2]. These effects originate from 
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an anisotropy in momentum space and lead to effective fields which are 
different from the ground state [3]. 

In this contribution we discuss the origin of non-equilibrium features, 
their influence on the nuclear matter EOS and the consequences when draw- 
ing conclusions on the nuclear EOS from comparisons with experiments in 
terms of collective flow effects. 

2. The nuclear EOS: from nuclear matter to heavy ion collisions 

The theoretical concept of the description of nuclear matter (NM) is given 
in an effective relativistic quantum field theory, Quantumhadrodynamics 
(QHD) [4]. It is formulated in terms of baryon and meson fields as the im- 
portant non-perturbative degrees of freedom of QCD. A relativistic treat- 
ment is attractive since it yields a natural description of characteristic fea- 
tures of NM, such as the appereance of strong attractive scalar and repulsive 
vector fields which naturally account for the saturation mechanism and for 
the strong spin-orbit potential. Relativity also leads to effectively energy 
dependent fields, which are crucial for nucleon-nucleus reactions and also 
heavy ion collisions. The effective Lagrangian of QHD involving baryon and 
meson fields is given by (for brevity, only scalar a and vector lo^ meson fields 
are considered here) 

£qhd = £b + -Cm + Ant 

Cm = \ (dpoWo - mla 2 ) - -F^ - ml w a u a 

C mt = iy^a - Tji la ^u a , (1) 

with F^ v = 8^uj u — d u uj fM the field tensor. The mesons a and vector uj^ are 
coupled in a minimal way to the nucleons via effective couplings T a and 
respectively. 

Different approximations have been used for the solution of eq. (1) in 
nuclear matter. The most popular one is the Hartree or Relativistic Mean 
Field (RMF) approach where the meson fields are treated classically [4]. 
A more realistic treatment is the Dirac-Brueckner-Hartree-Fock (DBHF) 
theory, in which exchange terms and higher order correlations within a 
T-Matrix or ladder approximation are taken into account [5]. The success 
of DBHF theory was in the unified description of nucleon-nucleon (NN) 
scattering and saturation properties of nuclear matter which was not pos- 
sible in non-relativistic approaches, except by including 3-body forces. The 
nucleonic mean field potential is given in terms of self energies by 



V DBHb (p,p F ) = Z s (p,p F ) - ^(p, PF ) 



(2) 
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Figure 1. Nuclear matter EOS for different models: (solid) DBHF from [8], (dashed) 
DBHF from [9], (long-dashed) RMF with density dependent couplings [10] and (dotted) 
Walecka model (NL2 parametrization) [11] 



with scalar and vector components S s , £ M depending on density p(pf) and 
for energies higher that the Fermi energy on momentum p. For NM the 
DBHF theory is parameter free, since it is based on a model for the bare 
NN interaction given by boson exchange potentials [6]. This is in contrast 
to phenomenological approaches, where several parameters are adjusted to 
NM saturation properties. 

The application of DBHF theory to finite nuclei has been formulated 
in a Density Dependent Hadronic (DDH) field theory [7], where the self 
energies are parametrized in Hartree form by 

S s = T s (p,p F )p s (p F ) , T lu = T (p,p F )j u , (3) 

with density and momentum dependent vertex functions r Sj o that replace 
the constant values of the RMF approach (r S) o — > T^/m^). They effec- 
tively contain exchange and correlation effects of the DBHF theory. Fig. 1 
compares the density dependence of the ground state EOS obtained from 
different models, which show a similar density behavior around saturation, 
but signifigant differences at high densities. The density dependence of the 
EOS at supra-normal densities can be regarded as an extrapolation which 
is tested in heavy ion collisions, where high compressions of the matter are 
reached. 

However, the determination of the nuclear matter EOS from heavy ion 
collisions is only indirect, because it contains fields corresponding to non- 
equilibrium colliding matter. Thus one does not, in fact, see the equilibrated 
nuclear matter EOS directly in heavy ion collisions. 

Equilibrated nuclear matter is characterized by an isotropic spherical 
momentum distribution. This is not the case in the dynamical situations 
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Figure 2. Time evolution of the central pressures (left figure) (longitudinal Pl and 
transversal Pr) and the central density (right figure) for a central An + Au heavy ion 
collision at Ei a b = 0.4 A.GeV 

of heavy ion collisions. This is seen in Fig. 2 in terms of the transversal 
and longitudinal components of the local central pressure, obtained in cal- 
culations with a transport equation, as described in the next section. The 
anisotropy in the pressure components indicates that the matter is locally 
not equilibrated, except at the late stage of the procces (t > 30 fm/c). 
Indeed, the relaxation times are large and comparable with the compres- 
sion phase of the collision (see time evolution of the central density in 2). 
Therefore, it is important to first study possible influences of phase space 
anisotropies on the nuclear matter EOS before making comparisons with 
experiments. 

The investigation of non-equilibrium effects on the level of the effective 
mean fields has been done by considering two models: a Local Density and 
a Colliding Nuclear Matter approximations, denoted as LDA and CNM re- 
spectively. In the LDA the EOS and the corresponding fields are those of 
equilibrated matter decribed by an isotropic momentum distribution. The 
mean field T,(pp,p) = T ls (pp,p) — 7 M S^(pi?,p) depends on a density p(Ep) 
and for energies above the Fermi energy Ep on momentum p relative to the 
rest system of nuclear matter. The density dependence has been taken from 
phenomenological (Walecka, DDH) or microscopic (DBHF) models. The 
momentum dependence is also naturally provided from the DBHF theory, 
whereas in the DDH approach it has been adjusted to the energy depen- 
dence of the Schrodinger equivalent optical potential in nucleon-nucleus 
scattering [12]. The fields of the non-linear Walecka model (NL2) do not 
depend on energy. 

In the CNM model the anisotropy of momentum space is modeled by 
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Figure 3. Shematic representation of the momentum space in the LDA and CNM models. 



that of two covariant Fermi-ellipsoids, i.e. two interpenetrating currents 
of nuclear matter (see Fig. 3) fcNM(pF 1 ,PF 2 ,v re l) = /i(PfJ + h{PF 2 ) + 
5f (pFi 2) v rel) where the last term takes Pauli effects into account. This 
representation of phase space anisotropies is motivated from microscopic 
calculations of heavy ion collisions [13]. The mean field is then given (e.g. 
for the scalar part) approximately by 

Yp NM (p, p p] Vrel ) = S s i (A^p, PFi)+S s2 (A2 1 P, PF 2 )+SZ(A^P, p Fl 2 , Vrel) 

(4) 

i.e. as a superposition of self energies of the nucleon with respect to the two 
currents taking into account the Lorentz transformations Ai 2 for the mo- 
menta and the Pauli correction in the overlap region (<5£). This ensures the 

correct limit to the LDA of one Fermi-sphere, i.e. CNM Vr ^L^° LDA. As 
a further approximation an average of the p-dependent quantities S^^^f,, 
over the CNM configuration leads to fields depending only on the configura- 
tion parameters, i.e. the two Fermi momenta pf 1/2 and the relative velocity 

Vrel [14]. 

The energy-momentum tensor T^ u is given in terms of the CNM fields 
by [3, 14] 



T au = < p *» p *v/ E * >CNM _ < 



CNMu 



/E* >CNM 



2 y 



< E,m7F >cnm - < P * X ^ NM /E* > C NM 



,(5) 



where < • • • >cnm means the average over the CNM momentum space. 
The energy per nucleon E CNM (ptot,v re i) = T 00 / p to t — M depends on the 
total invariant density p to t = Vltotajtot an d the relative velocity v re i- How- 
ever, the total energy E CNM contains contributions from the energy of 
the relative motion of the two currents. A meaningful discussion of non- 
equilibrium effects with respect to the ground state EOS should be based 
on the binding energy and thus it is reasonable to subtract contributions 
from the relative motion of the two currents E re i- E re i can be relativistically 
defined as the difference of the kinetic energy of the CNM system and the 
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Figure 4- Effective binding energy per particle fec/VM as function of the total density 
Ptot (normalized to the saturation density p sa t) at different relative velocities v re i (solid: 
Vrei = 0, dotted: v re i — 0.2, dashed: v re i — 0.4 and long-dashed: v re i = 0.6) for the DBHF 
theory from Refs. [8] (DB-HM) and [9] (DB-F) and for the DDH approach [10]. 



corresponding LDA configuration at twice the subsystem density 2p 



E re l(Ptot,Vrel) =< \f\ 



< \/p* 2 + m* 2 - 



■m >cnm 



< \/p* 2 + m* 2 - 



-m > v 



(6) 

By definition E re \ contains the kinetic energy arising from the separation 
of the Fermi-spheres in momentum space. It accounts for the interaction 
between the two currents by the presence of the effective mass in eq. (6). 
One should note that eq. (6) is the natural definition of a relativistic kinetic 
energy of the relative motion of two interacting currents and contains, in 
particular, the correct non-relativistic limit [3]. 

Now one is able to construct an effective bound energy of colliding 
nuclear matter E h ^f M {p toU v Td ) 



EcN d M(Ptot,V re l) 



T oo _ R 



rel 



M 



Ptot 



(7) 



which relates colliding to ground state nuclear matter. EQ^f M is displayed 
in Fig. 4 for different models. The case v re [ = is the correct limit for 
LDA. With increasing v re i we observe a softening of the effective EOS in 
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terms of an increase of the binding energy at p to t ~ Psat and a decrease 
of E^j^f M at high densities (ptot » Psat) with respect to the ground state 
EOS (solid curves in Fig. 4). The softening of the effective EOS occurs in 
all the models investigated and can be explained by (a) a decrease of the 
fields S^o^.. with density and momentum, see eq. (4), (b) a reduction of 
the Fermi pressure (not shown here) in transverse direction and (c) an en- 
larged scalar attraction with increasing relative velocity [3]. It is important 
that the differences between the non-equilibrium and the ground state EOS 
are of the same magnitude as the differences of the ground state EOS for 
different models, see e.g. Fig. 1. Thus we conclude that one should take 
non-equilibrium effects into account when determining the EOS in heavy 
ion collisions. In the next section we investigate the non-equilibrium effects 
in transport calculations of heavy ion collisions in terms of collective flow 
observables. 



3. Heavy ion collisions 

Heavy ion collisions have been described by covariant transport equations 
of a Boltzmann type (details can be found in refs. [11, 15] ): 

(m*d>* - p w d£pt) &l - (m*d£m* - p*"d£p* u ) d%] f(x, p) 

f(x,P 3 )f(x,P4) (1 - f(x,p)) (1 - f(x,p 2 )) ~ 

/(x,p)/(x,p 2 )(l-/(x,p3))(l-/(x,p 4 )) ] , (8) 

with effective quantities p* M = p^ — T,^ and m* = M — S s . It describes the 
evolution of the phase space density f(x,p) under the influence of a mean 
field in the left hand side of eq. (8) and on the right hand side of 2-body 
collisions determined by energy and isospin dependent NN-cross sections 
<tnn ~ W (pp2\P3P4) for the procces p + P2 — ► Ps +Pa- The Pauli principle 
is taken into account for the final states by the terms (1 — /). The NN-cross 
section is taken from empirical NN-scattering data. Also inelastic channels 
are taken into account as the production of A and ./V* resonances with their 
decay in 1- and 2-pion channels [2, 16]. The EOS enters into the theoretical 
descriptions of heavy ion reactions via the self energies. 

We have applied the approximations discussed in the previous section, 
the LDA and CNM approaches, in transport calculations using eq. (8). In 
the LDA model the ground state EOS enters directly into the transport 
equation neglecting non-equilibrium effects, whereas in the CNM model 
the effective EOS including the momentum space anisotropies is consid- 
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Figure 5. Energy dependence of the elliptic flow for peripheral Au + Au collisions at 
energies from the Fermi energy to AGS energies. Calculations with the DBHF model [9] 
are shown in the LDA and CNM approximations and compared to data [17]. 



ered in eq. (8). In the latter case the CNM parameters, i.e. the invariants 
Pi-, P2-, v re i, are determined directly from the phase space density f(x,p) 
at each position and in each time step of the simulation. 

We discuss the results of transport simulations in terms of collective 
flows, which have been found to be sensitive to the density and, partic- 
ularly, momentum dependence of the mean field [18]. The collective flow 
is characterized by a Fourier series of the azimuthal distribution of nucle- 
ons for given rapidity y and total transverse momentum p t , N(cj),y,pt) = 
Nq(1 + v\coscj) + V2Cos2(ft • • •), where v\ and V2 are called transverse and el- 
liptic flow, respectively. As an example, Fig. 5 shows the energy dependence 
of the elliptic flow at mid-rapidity (|Ay| < 0.15) for peripheral Au + Au 
collisions. With the cut in the rapidity y = E ^ J Z one selects particles from 
the hot and compressed matter in the central region of the reaction, from 
which one intends to extract information on the nuclear EOS. Thus the el- 
liptic flow describes the dynamics of the compressed fireball perpendicular 
to the beam direction. It probes the high density behavior of the nucleonic 
fields, i.e. their stiffness, which gives rise to a squeeze-out of the compressed 
matter perpendicular to the beam direction. An enhancement of the elliptic 
flow is correlated to a stiffer EOS. The calculations in Fig. 5 were performed 
with one of the models of the previous section, i.e. the DBHF results of the 
Tubingen group [9] treated in the LDA and in the CNM approximations. 
The two different treatments which are based on identical nuclear forces for 
ground state matter yield significantly different results for the elliptic flow. 
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Figure 6. Transverse momentum (pf^ = \J 'pi + p\ IpY°^ ) dependence of the collective 
flow projected in (v\) and out of (vi) the reaction plane. (Left) Comparisons for Ru + Ru 
and (right) for Au + Au collisions at 0.4 A.GeV beam energy. The figure on the right 
shows the quantity vi for protons (Z — 1) and fragments (Z = 2). The data (symbols) 
are from ref. [17] and the theoretical calculations were performed within the DBHF theory 
of refs. [9, 8] in the CNM model. 



Since this observable is particularly sensitive to the EOS at high densities 
[19, 18], these effects are most pronounced here. In the LDA approach one 
would have excluded the underlying EOS from the comparison to data as 
too stiff. The more consistent treatment of momentum space anisotropies 
on the level of the effective interaction leads to the net softening of the 
effective EOS discussed above and restores the agreement with experimen- 
tal data. Similar effects are observed for other models, e.g. those of the 
previous section, and other components of collective flow not shown here 
[2, 16]. 

Other examples of collective flow effects are shown in Fig. 6 with respect 
to the transverse momentum dependence pf^ . High p^' -values correspond 
to highly energetic particles originating from the compressed matter in an 
early stage of the reaction and thus directly carry information about the 
high density behavior of the EOS. The DBT calculations with a relatively 
soft EOS also yield a smoother p^-dependence of the collective flows closer 
to the data than the DBHM calculations. This effect, which is observed 
also for other energies, is consistent with Figs. 4 and 5 indicating that a 
moderate dependence of the nuclear matter EOS at high densities seems 
more realistic (similar effects have been found with the phenomenological 
DDH model). 
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4. Conclusions 

The early and high density phases of relativistic heavy ion collisions are 
to a large extent governed by highly anisotropic phase space configura- 
tions far from local equilibrium, which can be approximated by counter- 
streaming or colliding nuclear matter configurations. In this contribution we 
discussed the implications for the effective EOS which occurs in such non- 
equilibrium configurations and found them to be important. The anisotropy 
in momentum space leads to a softening of the non-equilibrium EOS with 
essential deviations compared to the ground state EOS. CNM calculations 
with different effective interactions for the mean field have shown the gen- 
eral features of this effect. In calculations of energetic heavy ion collisions, 
where the mean field were considered in the local density (LDA) and in 
the more realistic CNM approximations, the effective softening of the non- 
equilibrium EOS was observable in terms of collective flow effects. The 
differences in the flow signals are found to be comparable to those of the 
nuclear matter EOS for different models. We conclude that for a reliable 
determination of the nuclear matter equation of state from heavy ion col- 
lisions the non-equilibrium features of the phase space must be taken into 
account on the level of the effective fields. The results seem to be consistent 
with a moderate dependence on momentum of the nuclear matter EOS at 
high densities. With the successful applications of microscopic models to 
nuclear matter, finite nuclei and heavy ion collisions one expects to move 
towards a unification of the description of very different nuclear systems 
and to a determination of the EOS of nuclear matter. 
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